arXiv:1509.01794vl [physics.atom-ph] 6 Sep 2015 


Theoretical studies of the long lifetimes of the 6 d ^i?3/2,5/2 states in Pr: Implications for 

parity nonconservation measurements 

B. K. Sahoo * 

Theoretical Physics Division, Physical Research Laboratory, Ahmedabad-380009, India 

B. P. Das 

International Education and Research Center of Science and Department of Physics, 

Tokyo Institute of Technology, 2-1-2-1-H86 Ookayama Meguro-ku, Tokyo 152-8550, Japan 
(Dated: Recieved date; Accepted date) 

The lifetimes of the 6d ^D^j^ and 6d ^ 1 ) 5/2 states in Fr are estimated to be 540(10) ns and 1704(32) ns 
respectively. They are determined by calculating the radiative transition amplitudes of the allowed 
electric dipole (El) and the forbidden electric quadrupole (E2) and magnetic dipole (Ml) channels 
using the second order many-body perturbation theory (MBPT(2)) and the coupled-cluster (CC) 
method at different levels of approximation in the relativistic framework. These long lifetimes and 
the large electric dipole parity non conserving amplitudes of 7s ^S' 1/2 6d ^D^/ 2 ,s /2 transitions 
strongly favour Fr as a leading candidate for the measurement of parity nonconservation arising 
from the neutral current weak interaction and the nuclear anapole moment. 

PACS numbers: 


Francium (Fr) is considered to be a promising candi¬ 
date for the measurements of the electric dipole moment 
(EDM) arising due to the violations of parity and time 
reversal symmetries Hi , parity nonconservation (PNC) 
effects due to the neutral weak interaction ii and the 
nuclear anapole moment Si as it is the heaviest alkali 
metal atom. The focus of all the ongoing Fr PNC ex¬ 
periments is the 7 s ^S'1/2 —>■ 8s ^51/2 transition Sii; 
which is largely inspired by the Cs PNC experiment g. 
However, relativistic many-body calculations show that 
the PNC amplitudes in the 7 s ^S'1/2 ^ Sd ^^3/2,5/2 tran¬ 
sition amplitudes in Fr are about three times larger than 
that of the 7 s ^S'1/2 ^ 8s ^S'1/2 transition Si- The 
S — D PNC transitions in the singly ionized Ba, Ra and 
Yb have been the subject of theoretical investigations 
id and the principle of their measurements has been 
discussed [ 1 ^ |l^ . It has also been highlighted that the 
PNC measurements for the S — D5/2 transitions in these 
ions would provide unambiguous signatures of the exis¬ 
tence of the nuclear anapole moment (NAM) Bin, 
which is still an open question [l^ . Apart from exhibit¬ 
ing large PNC effects, another important aspect of these 
transitions is that the excited D states in theses ions 
are metastable stables, and they provide long interro¬ 
gation times which enhances the precision of the mea¬ 
surement of the small PNC effects [l^. In this Rapid 
Communication, we present the results of our theoreti¬ 
cal studies on the lifetimes the 6d ^^3/2 and 6d "^0^/2 
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states in Fr, which were undertaken to assess the feasi¬ 
bility of the measurement of PNC in this atom using the 
7 s ^5 'i/ 2 ^ 6d ^^3/2,5/2 transitions. 

An electron from the 6d ^£>3/2 state can decay to its 
low-lying 7 p ‘^Pi/2 and Ip ^^3/2 states by the electric 
dipole (El) and forbidden magnetic octupole (M 3 ) tran¬ 
sitions and to the ground state by the forbidden magnetic 
dipole (Ml) and electric quadrupole (E 2 ) transitions. We 
neglect contributions due to the M 3 transition as the cor¬ 
responding transition probability is very weak owing to 
its inversely proportional to seventh power of transition 
wavelength. Similarly, an electron from the 6 d ^1)5/2 
state can decay to its fine-structure partner 6d ^D3/2 
state via both the Ml and E 2 transitions while to the 
low-lying Ip ^^3/2 state by the El transition and to the 
ground state by the E 2 transition. In this case too, we 
have omitted contributions due to the M 3 transition. The 
transition probabilities due to the above El, E 2 and Ml 
channels for a transition, say, | 5 'i) ^ |'I'/) are given by 

3B1 _ 2.0261 X 10 ® ^£3 

~ \3 „ 

\f9^ 

aE 2 _ 1-1195 X 10 (2_E2 

\f9^ 

and 

Ml _ 2.6971 X 10 ^Mi 

^if ~ \3 > 

\f9i 

where the quantity =| ('I'i||0||T'/) ^ is known as the 
line strength for the corresponding reduced matrix ele¬ 
ment I {T'i|| 0 || 4 'y) I of a transition operator O. These 
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quantities are later given in atomic unit (a.u.). In the 
above expressions, gi = 2 + 1 is the degeneracy factor 

of the state Idti) with the angular momentum Ji and the 
transition wavelength (A^/) is used in nm which when 
substituted, the transition probabilities {A^^s) are ob¬ 
tained in The lifetime (r) of the atomic state |d>i) is 
determined by taking the reciprocal of the total emission 
transition probabilities involving all the possible sponta¬ 
neous transition channels (in s). i.e. 


Ti = 


1 




(4) 


where the summations over O and / correspond to all 
probable decay channels and all the lower states respec¬ 
tively. We have attempted to obtain accurate results 
for the transition probabilities, hence the lifetimes of the 
atomic states by performing relativistic many-body cal¬ 
culations of the line strengths and using wavelengths that 
are determined from the experimental transition energies 
given in the National Institute of Science and Technology 
(NIST) database [l^ . 

To investigate the role of the electron correlation effects 
in the evaluation of the radiative transition amplitudes, 
we employ the second order many-body perturbation the¬ 
ory (MBPT(2)) and the coupled-cluster (CC) method in 
the relativistic framework. Further, we take different lev¬ 
els of approximation in the CC method to see the conver¬ 
gence in the results. We give below a brief description of 
these methods using Bloch’s prescription (2^, in which 
atomic wave function of state IT^) is expressed as 


(5) 

where and |$„) are known as the wave operator and 
reference state respectively. The ground and the consid¬ 
ered excited 6D states of Fr have the electronic configura¬ 
tions as [6p®]7s ^Si /2 and [6p^]6d ^ 113 / 2 , 5/2 respectively. 
To reduce the computational effort, we construct these 
states by creating a common reference state function |$c) 
with the [6p®] configuration using the Dirac-Hartree-Fock 
(DBF) method. In this approach the atomic Hamiltonian 
(H) in the Dirac-Coulomb interaction approximation is 
divided as DHF Hamiltonian Hq and residual Coulomb 
interaction I 4 . For the calculation of the exact states 
with a valence orbital, we define new working reference 
states as |$„) = Here aj^ appends an electron 

from the respective valence orbital denoted by an index 
n. As a consequence can now be divided as 


In these expressions, the superscripts imply number of 
Vr considered in the calculations and represents order of 
perturbation; e.g. MBPT(2) method bears terms up to 
two Vr {k = 2). 

Using the generalized Bloch’s equation, fcth order am¬ 
plitudes for the Xc and Xn operators are obtained by [ 2 ^ 

= QVr{lVx^t^^)P ( 8 ) 

and 

k-l 

[xi"\ = QVr{l + + X^t^^)P - E 

m—1 

xPK(I + xi™-')+xi’""'^)i", (9) 

where the projection operators P = |d)c)(d)c| and Q = 1 — 
P describe the model space and the orthogonal space of 
the Hamiltonian Hq respectively. Note that here Xc°^ = 0 
and Xn^ = 0- Using these amplitudes, the energy of the 
state is evaluated by using an effective Hamiltonian 

= PHn^P ( 10 ) 

Using the CC ansatz, the above expressions can be put 
together to construct a wave operator to infinite order as 

|«-„) = =e^{l + 5„}|$„), (II) 

such that Xc = — 1 and Xn = — I. Here the 

T and Sn are the CC excitation operators that excite 
electrons from the core and core along with the valence 
orbital respectively. In this work, we have only accounted 
the singly and doubly excited states denoting the CC 
operators by subscripts 1 and 2 respectively as 

T = Ti+T 2 and + ^ 2 ™. (12) 

This is referred to as the CCSD method in the literature. 
When only the linear terms are retain in Eq. (TO with 
singles and doubles approximation, it is referred to as 
the LCCSD method. The amplitudes of the T and Sn 
operators are determined using the expressions 

HnXcP = QHnP (13) 


and 


HnXuP = QHN{l + Xc)P-XnH‘^J^. (14) 


Un = l+Xc + Xn, (6) 

where Xc and Xn are responsible for carrying out excita¬ 
tions (generating configuration state functions) from |$c) 
and |$„), respectively, due to W- In a perturbative series 
expansion, we have 

Xc = EAi"^ and Xn=EAi'=^- (7) 

k k 


where we have defined normal order Hamiltonian Hn = 
H — PHP and the effective Hamiltonian to evaluate ion¬ 
ization potential (IP) of an electron from the valence 
orbital n of the respective state is given by = 

PHn{1 + Xc + Xn)P- We have also included contribu¬ 
tions from important triples excitations perturbatively 
from Xc^^ and Xn^ in the construction of and the 

approach is referred to as the CCSD(T) method. 
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TABLE I: Magnitudes of the reduced matrix elements {Ji\\0\\Jj) of transition operators (Os) given in a.u. from different 
many-body methods. Uncertainties from the finite size basis set, non-inclusion of the Breit interaction and due to QED effects 
are quoted using the MBPT(2) method. Wavelengths (Ai/) from the NIST database 0 are quoted in nm for the respective 
transitions. 

Ji -5- Jf Xif DHF MBPT(2) LCCSD CCSD^^' CCSD'-^) CCSD(°°) CCSD(T) Uncertainties 


Basis Breit QED 


El matrix elements 


6d ^D^J2 7p ^Pi/2 

2504.7 

9.22 

7.73 

6.81 

7.46 

7.47 

7.47 

7.43 

0.05 

-0.01 

0.001 

6d ^D^/2 7p ^P3/2 

4336.7 

4.28 

3.57 

3.12 

3.44 

3.45 

3.45 

3.42 

0.02 

-0.01 

~ 0.0 

6d ^D^/2 7p ^^3/2 

3991.0 

12.80 

10.83 

9.68 

10.54 

10.55 

10.55 

10.51 

0.07 

-0.02 

0.001 


E2 matrix elements 


6d ^Sij2 

616.5 

43.10 

33.74 

31.39 

34.02 

34.06 

34.06 

33.78 

0.20 

-0.03 

-0.02 

6d ^D^^2 ^'S’ 1/2 

608.7 

52.74 

41.69 

39.31 

42.24 

42.27 

42.27 

41.96 

0.24 

-0.04 

-0.02 

6d ^D^I2 Qd ^d^j2 

50057.6 

47.70 

32.14 

27.18 

32.01 

32.03 

32.03 

31.49 

0.55 

-0.16 

0.01 





Ml matrix elements 






6d ^D^f 2 7s ‘^S\j 2 

616.5 

~ 0.0 

0.001 

0.002 

0.001 

0.001 

0.001 

0.001 

- 0.0 

- 0.0 

~ 0.0 

6 (^ ^D^/2 ^t/ 3/2 

50057.6 

1.549 

1.550 

1.548 

1.547 

1.552 

1.552 

1.552 

~ 0.0 

0.001 

~ 0.0 


After obtaining amplitudes using the above equations, 
the transition matrix element of an operator O between 
the states Idfi) and Idt/) is evaluated using the expression 


the CCSD^^) method considering terms up to /c ^ oo 
in a self-consistent procedure as described above. Our 
final CCSD results correspond to the CCSD(°°) method. 
The same procedure is also adopted for the CCSD(T) 
method. The contribution from the normalizations of 
the wave functions (Cnorm) is estimated explicitly using 
the expression 


This gives rise to a finite number of terms for the 
MBPT(2) and LCCSD methods, but it involves two non¬ 
terminating series in the numerator and denominator 
which are and respectively in the CCSD 

and CCSD(T) methods. In order to evaluate all the sig¬ 
nificant contributions from these series, we have used the 
Wick’s generalized theorem to divide these terms into 
the effective one-body, two-body and three-body terms. 
The effective one-body terms are the dominant ones, they 
are computed first considering the CC terms with the ap¬ 
proximations Oe'^ ~ 0+0T+T'^0+^0T'^ + ^T^'^0+ 
T'lQT and ~ T^T + -b Then, they 

are stored and contracted with the T 2 and operators 
avoiding repetitions of the diagrams in a self-consistent 
procedure to account for the higher order one-body terms 
from the non-terminating series. They are again stored 
as an intermediate form for the further contraction with 
the Sn and operators. Similarly the effective two- 
body and three-body terms are computed after contract¬ 
ing with the above effective one-body terms with the T 2 
and T 2 operators, but they are computed directly con¬ 
tracting with the Sn and S'! operators. Thus, these effec¬ 
tive two-body and three-body terms also have contribu¬ 
tions from the non-terminating series. To see the conver¬ 
gence of the results with the series expansion, we present 
contributions with k numbers of T and/or operators 
from these non-terminating series, which we refer to as 


C 


norm 






(16) 


In Table HI we give the radiative transition matrix el¬ 
ements for all the considered channels from the DHF, 
MBPT(2), LCCSD, CCSD and CCSD(T) methods to an¬ 
alyze the propagation of the correlation effects through 
various levels of approximations in the many-body theo¬ 
ries and the experimental values of the transition wave¬ 
lengths from the NIST database 0 that we have used 
later. We also give contributions from the CCSD method 
by truncating the non-linear terms with fc = 2, fc = 4 and 
from a self-consistent {k = 00 ) calculation. For fc = 2, the 
expression evaluating the property given by Eq. (113 has 
the same number of terms as does the LCCSD method. 
Therefore, differences in the results from the LCCSD 
and CCSD^^) methods imply the correlation contribu¬ 
tions arising through the non-linear terms in the wave 
function determining equations of the CCSD method and 
are found to be quite large. Often, these contributions 
are neglected in the calculations as they require pro¬ 
hibitively large computational resources for their eval¬ 
uation. Again, we observe from the trends that the cor¬ 
relation effects at the MBPT(2) method are large, and 
that there are strong cancellations in the LCCSD approx¬ 
imation and the results almost converge for k = A when 
non-linear terms are included in the CCSD method. The 
















4 


TABLE II: Contributions to the reduced matrix elements {Ji\\0\\Jf) from various terms of the CCSD method (in a.u.). 
Differences between these values from the CCSD results are quoted in Table U correspond to those non-linear terms that are 
not mentioned explicitly here. 


Ji —y Jf 


0 

OTi 

t/o 

05i/ 

5^0 

5t,05i/ 

052/ 

slo 

51.052/ 

5^051/ 

C-norm 






El matrix elements 






6d ^D^/2 

—>■ 7p ^Pl /2 

9.22 

- 0.0 

0.019 

-0.437 

-0.877 

0.246 

-0.208 

-0.248 

- 0.020 

- 0.020 

-0.324 

6 d ^D^/2 

—>■ 7p ^£ 3/2 

4.28 

~ 0.0 

0.002 

-0.132 

-0.522 

0.089 

-0.081 

-0.097 

-0.005 

-0.007 

-0.139 

6d ^£> 5/2 

7p ^£ 3/2 

12.80 

~ 0.0 

0.006 

-0.381 

-1.372 

0.238 

-0.233 

-0.292 

-0.008 

- 0.020 

-0.342 






E2 matrix elements 






6 d ^0^/2 

^ 7s 251/2 

43.10 

~ 0.0 

0.013 

-6.306 

-2.683 

1.547 

-0.139 

-0.111 

-0.027 

-0.186 

-1.618 

6d ^D^/2 

^ 7s 251/2 

52.74 

~ 0.0 

0.015 

-7.644 

-2.871 

1.666 

-0.125 

-0.148 

-0.017 

0.004 

-1.6 

6d ^0^/2 

—>■ 6d 2^3/2 

47.70 

- 0.0 

~ 0.0 

-8.994 

-8.116 

3.264 

-0.079 

-0.010 

-0.012 

-0.021 

-1.947 






Ml matrix elements 






6d ^0^/2 

^ 7s 251/2 

~ 0.0 

- 0.0 

~ 0.0 

- 0.0 

- 0.0 

0.0004 

0.0005 

0.0001 

- 0.0 

- 0.0 

- 0.0 

6d ^D^/2 

6d 2(13/2 

1.549 

~ 0.0 

~ 0.0 

0.003 

-0.002 

0.078 

0.003 

- 0.0 

~ 0.0 

~ 0.0 

-0.093 


TABLE III: Transition probabilities (A^) due to different 
transition decay channels (Os) in and their branching 
ratios from the 6d ^D ^/2 and 6d ^D ^/2 states of Fr. Uncer¬ 
tainties are quoted within the parentheses. 


Ji — y Jf 

0 

Aif 

Branching ratio 

6d ^D^/2 7p ^Pi/2 

El 

1779511(33688) 

0.96 

Qd ^0^/2 7p ^£3/2 

El 

72637(1280) 

0.04 

6d ^D^/2 —>■ 7s ^5i/2 

E2 

35.96(60) 

- 0.0 

6d ^D^/2 7s ^5i/2 

Ml 

~ 0.0 

- 0.0 

6d ^D^/2 7p ^£ 3/2 

El 

586776(11219) 

- 1.0 

Qd ^'S’l/2 

E2 

39.33(19) 

~ 0.0 

Qd ^D^j2 Qd ^d^j2 

E2 

- 0.0 

- 0.0 

Qd ^D^^2 ^^3/2 

Ml 

- 0.0 

~ 0.0 


discrepancy in the results of the CCSD*^^) and CCSD(°“) 
methods is beyond second significant digit implying that 
the results have converged within the precision of our 
interest. The valence triple excitations seem to change 
the results slightly. We also give uncertainties associated 
with these results by estimating contributions due to the 
finite size of our basis set, neglected contributions from 
the Breit interaction and corrections from the quantum 
electrodynamics (QED). These estimates are carried out 
using the MBPT(2) method which gives the largest cor¬ 
relation effects. 

After analyzing the trends in the correlation effects at 
different levels of approximation, we now focus on the 
contributions from different terms in the CCSD method. 
We present these results in Table HIl along with the contri¬ 
butions from Cnorm- Contributions from the correspond¬ 
ing radiative operator O are the DHF results, OTi and its 
complex conjugate terms give core-valence correlations, 
OSif and Slfi give the pair correlation effects involv¬ 
ing the valence orbitals, 052 / and Sl^O give the core¬ 


polarization correlation effects involving the valence or¬ 
bitals, etc. Contributions from the other non-linear terms 
such as representing the core pair correlation effects com¬ 
ing through the T 2 OT 2 term are not given explicitly in 
the above table, however their contributions can be ob¬ 
tained by taking the differences of the contributions given 
in Table inland the final CCSD results given in TableH] As 
can be seen from Table m core-valence correlations are 
small and the largest correlation effects come through the 
pair correlation effects in the El and E2 matrix element 
calculations. Nevertheless, the core-polarization effects 
are also very significant and are the dominant ones in 
the calculations of the Ml matrix elements. We also find 
contributions from Cnorm to be fairly large. 

We now use the matrix elements from the CCSD(T) 
method and the experimental wavelengths mentioned in 
Table U to evaluate the transition probabilities due to 
different radiative decay channels from the 6d ^^ 3/2 and 
6d ^D ^/2 states. These values are given in Table Hill along 
with their branching ratios for the individual transitions. 
As can be seen from this table that branching ratios are 
dominated by the El transitions and they are entirely 
responsible for determining the lifetimes of the 6d ^ 0^/2 
and 6d '^D^i 2 states. Using these transition probabilities, 
we estimate the lifetimes of the 6d and 6d ^£> 5/2 

states as Ted 3/2 = 540(10) ns and TQd 5/2 = 1704(32) 
ns, respectively. These values are very large compared 
to the other low-lying excited 7p ^£ 1 / 2 , 372 , 8 s ^ 5 i/ 2 , 
7p ^£’i/ 2 . 3 / 2 j 7d ^D 3 / 2 , 5/2 and 9s ^ 5 i /2 states of Fr which 
are measured till date as T'j'pi /2 = 29.45(11) ns, T’jp 3 i 2 = 
21.02(15) ns, r 8 ,i /2 = 53.30(44) ns, T^a 3/2 = 73.6(3) 
ns and T’jii ^/2 = 67.7(2.9) ns, Tspi /2 = 149.3(3.5) ns, 
T‘ 8 p 3/2 = 83.5(1.5) ns and T 9 S 1/2 = 107.53(90) ns re¬ 
spectively [llj. Given the large PNC amplitudes in the 
7s ^ 5 i /2 6 d ^ 03 / 2 , 5/2 transitions 3 compared to the 
7s ^ 5 i 72 ^ 8s ^ 5 i /2 transition in Fr M and with the cor¬ 
responding transitions in Ra'*' ion [13, HBl and the long 
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lifetimes of the excited 6d ^ 0^/2 and 6d ^D ^/2 states, 
make Fr a potentially attractive candidate for a PNC 
experiment. 

In summary, we have performed relativistic many- 
body calculations of the lifetimes of the 6d ^I? 3/2 and 
6d ^£> 5/2 states of Fr and they are found to be large. 
These results favour the measurement of PNC in the 
7s ‘^Si ^2 —^ bd ‘^D ^^2 and 7s ^5*3^2 —^ bd ‘^D^i 2 transi¬ 
tions of Fr, where calculations predict large PNC effects. 
For the evaluation of the lifetimes, we have calculated 
radiative transition matrix elements using the relativis¬ 


tic CC method. We have also investigated the roles of 
the electron correlation effects in the determination of 
these quantities systematically by approximating many- 
body methods at different levels and give contributions 
explicitly from various CCSD terms. 

This work was supported partly by INSA-JSPS un¬ 
der the project no. lA/INSA-JSPS Project/2013- 
2016/February 28,2013/4098. Computations were car¬ 
ried out using the Vikram-IOOTF HPC cluster at the 
Physical Research Laboratory, Ahmedabad, India. 
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